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ABSTRACT 

Zonally-averaged  autoeovariances  of  the  five-day  height  tendency  for 
lags  0,  5,  ,..,  25  days  are  computed,  at  25  mb  aro\ind  latitudes  55N,  65N, 
and  75N  during  the  period  5  October  1958  through  29  March  1959.  Using 
the  convolution  theorem  for  each  space-wave  number,  it  is  possible  to 
determine  the  zonally-averaged  spectral  density  of  the  height-tendency 
field  at  each  of  these  latitudes  in  terms  of  time  leaves  having  periods 
from  10  to  50  days.  The  results  indicate  that  all  but  about  10/.'.  of  the 
variance  of  the  height-change  field  is  accounted  for  by  space  waves  one 
and  two.  For  space  wave  number  two,  time  waves  number  3  and  4-  (per  50 
days )  give  a  significant  contribution  to  the  spectrum  at  the  90%   confi- 
dence level.  However,  for  space  wave  number  one,  time  wave  number  one 
gives  the  most  significant  contribution  to  the  spectrum.  It  is  inferred 
from  this,  that  cases  of  "explosive"  warming,  which  may  occur  as  early 
as  January  are  correlated  with  intensification  and/or  movement  associ- 
ated with  space  wave  number  one,  at  periods  of  approximately  100  days. 

The  author  wishes  to  express  his  sincere  appreciation  to  Professor 
Prank  L.  Martin,  Department  of  Meteorology,  IT.  3.  Naval  Postgraduate 
School,  for  his  suggestion  of  the  topic  and  his  encouragement  in  the 
preparation  of  this  study. 
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1 .  Introduction 

During  the  past  ten  years,  meteorologists  have  become  aware  of  the 
winter-time  phenomenon  known  as  "explosive"  warming,  which  occurs  in  the 
Arctic  and  sub-Arctic  stratosphere.  During  this  same  decade,  the  increased 
availability  of  higher-level  soundings  has  made  accessible  to  the  meteor- 
ologist sufficient  information  to  analyze  the  structure  and  causes  of 
these  unusual  occurrences. 

The  work  of  the  McC-ill  University  group,  primarily  under  F.  K.  Hare 
ll 960 1  and  B.  W.  Boville  n96o]  has  made  it  clear  that  there  is  a  succes- 
sion of  cold  troughs  and  warm  ridges  at  the  25-mb  level  which  eventually 
culminates  in  the  development  of  an  intense  warm  anticyclone  in  polar 
latitudes,  northern  Canada  is  a  preferred  location  for  the  occurrence  of 
this  phenomenon.  This  intense  anticyclogenesis  may  be  explained  |_Boville, 
196CJ  in  terms  of  baroclinic  instability  in  the  polar  night  westerlies. 
Actually  the  development  of  the  intense  anticyclone  represents  a  stage  in 
the  weakening  of  the  polar  night  westerlies  [Julian,  1 959J  which  are  re- 
placed in  spring  and  summer  by  the  easterly  wind  regime  that  finally  en- 
velops the  entire  summer  stratosphere.  Boville  and  bis  associates  [l96ll 
have  prepared  an  atlas  of  the  25-mb  circulation  for  the  six-month  period 
5  October  1958  through  29  March  1959.  In  these  studies,  Boville  et  al. 
list  the  Fourier  analyses  of  the  25-mb  contour  field  at  five-day  intervals 
for  latitudes  55N,  65N,  and  75N.  For  each  of  these  latitude  circles,  they 
describe  the  Fourier  analyses  in  the  form 

z  =  I  +2c„cos(nA-(l)n)       <u 

n=l 

where  the  values  of  C^,  and  wn   are  listed  in  the  atlas  through  space-wave 

number  nine.  Here  "2.  is  the  zonal  mean  contour  height  around  the  latitude 
circle  while  C„  and  (J)  are  the  amplitude  and  phase  angle,  respectively, 

1 


of  the  nth  harmonic  (n=1 ,....., 9 ) . 

Spectral  analysis  has  become  a  useful  diagnostic  tool  of  meteorolo- 
gists, primarily  for  the  purpose  of  detecting  preferred  periodicities  in 
the  eddies  or  perturbations  under  consideration.  The  number  and  variety 
of  such  studies  in  certain  areas  of  statistical  meteorology  are  too  great 
to  attempt  to  describe  them  here.  Suffice  it  to  say  that  the  basic  methods 
were  largely  developed  by  J.  W.  Tukey,  and  a  concise  summary  of  their  ap- 
plications to  Meteorology  can  be  found  in  the  excellent  textbook  by 
Panofsky  and  Brier  [l958]. 

Actually,  Boville  [i960]  has  published  the  spectral  analysis  of  the 
25-mb  contour  height  field  at  latitude  60N  for  the  month  of  October  1958. 
In  that  analysis,  one  of  the  conclusions  reached  was  that  most  of  the 
variance  in  the  contour  height  field  at  25  mb  is  explained  by  space-wave 
numbers  one  and  two.  A  similar  conclusion  is  reached  in  this  paper  in 
regard  to  the  variance  of  the  height  tendency. 

The  objective  of  the  present  paper  may  be  understood  from  the  local- 
time  differentiation  of  equation  (1 ) 


-^L   Pn  COS(n^"W+  Cn<j)n6iN(nA-(t)n)] 


(2) 


where   j?*  ss  <3  Z./c)t  —  ()£/  OL   •  From  equation  (2),  it  may  be 
seen  that  time-changes  in  the  perturbation  field  may  be  ascribed  to  terms 
involving  Cn,  the  local-time  change  in  the  harmonic  wave  amplitude,  and 
(J)n  ,  the  local-time  change  of  the  phase  angle.  These  two  types  of  terms 
in  equation  (2)  may  be  described,  as  the  intensification  and  movement  con- 
tributions, respectively,  for  the  nth  wave.  A  time-spectral  analysis  of 

d^/dt,   zonally  averaged  with  respect  to  longitude,  was  then  performed 
using  the  five-day  rates  of  change  C^  and  (pn .   (For  the  mathematical  de- 
velopment see  section  2.)  An  advantage  of  the  spectral  analysis  of  the 
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height-tendency  field  over  that  of  the  height  field  alone  is  that  the 
former  can  indicate  periodicities  due  to  intensification  and  to  movement 
of  the  harmonic  waves,  whereas  these  connections  are  not  so  clear  in  the 
latter  case. 


2.  Mathematical  development 

The  autocovariance  of  Z*  at  lag  t    is  given  [Blackman  and  Tukey, 
1958]  by  the  formula 

R(T)=^i/  ilt)  Z.*(t+T)dt  (3) 

J-T/Z 

where  Z*  is  the  zonal  height  perturbation,  and  includes  the  contribution 

from  all  space  waves  around  the  latitude  circle,  T  is  the  averaging  per- 
iod over  which  the  autocovariance  is  determined,  X   is  the  lag,  and  t  is 
an  interval  of  time,  also  measured  in  five-day  units.  Using  the  truncated 
time  series  available  for  this  study,  the  time  integration  may  be  replaced 
by  a  summation  over'ft-K  equally-spaced  time  intervals,  at  lag  ?=  KAt  . 
Then  the  integral  (3)  may  be  rewritten  in  the  form 

27T  ^H< 

R(K)=iJ  A  Iz*(t)  't(ux)]  M  o^K^vi ,     u) 

0        t=i 
where  R(k)  is  the  zonally-averaged  time  autocovariance.  In  equation  (4.), 

A  represents  the  longitude  at  any  point  within  a  latitude  circle.  Here,  K 
is  the  number  of  lags  in  terms  of  intervals  of  Al  =  5  days,  ^  is  the  num- 
ber of  five-day  height  changes  available  in  the  time  series,  and  M  is  the 
maximum  number  of  time  lags  used  in  this  study.  Obtaining  the  five-day 
perturbation  height  tendency  g^ft  +  T )     at  time  t+T   by  an  equation 
analogous  to  (2),  we  may  write, 

£*^)=H(lco5(nA-rt+&ijl  WnA-ij),,)}  (5) 

n-t 
where  Qr    and  (JJj,  are  analogous  to  Cn    and  $h  of  equation  (l).  Hence, 

from  equations  (2),  CO,  and  (5),  the  zonally-averaged  time  autocovariance 

for  space  wave  n  has  the  form 


[■ 


•  •  -i 


aQ^-^G^  I  SIN  ((U„  -(f),,)  J  .    (6) 

In  equation  (6)  the  subscript  n  indicates  that  RqOO  is  independent  of  all 
space  waves  other  than  the  nth.  In  the  derivation  of  (6),  considerable 
simplification  has  been  made  possible  as  a  result  of  orthogonality  prop- 
erties, that  is  by  deleting  terms  which  vanish  upon  integration  around  a 
complete  latitude  circle.  A  proof  of  equation  (6)  is  given  Appendix  I. 
Since  discrete  data  and  a  finite  time  series  have  been  used,'  the  maximum 
lag  MAX   limits  the  maximum  resolvable  period  to  be   |^ax  =  2.  liA  L  , 
which  may  be  called  the  fundamental  period.  In  this  study  the  maximum 
number  of  lags  has  been  taken  as  M  =  5  so  that  the  fundamental  period  will 
be  50  days.  If  K  is  allowed  to  range  over  the  values  0,  1,  2,  . .♦,  M, 
various  values  of  R^K)  can  be  generated  for  the  different  space-wave 
numbers.  Blackman  and  Tukey  [l958j  recommend,  when  using  a  truncated 
time  series,  that  the  maximum  number  of  lags  M  should  not  exceed  approx- 
imately 10$  of  the  length  of  the  record.  Since  the  total  available  number 
of  five-day  mean  values  of  ba/QX    was  35,  the  maximum  number  of  lags 
was  taken  to  be  five.  This  lag  number  was  considered  to  be  in  reasonable 
accord  with  the  recommendation  of  Blackman  and  Tukey,  and  yet  sufficiently 
large  to  give  resolution  for  the  long-period  waves  which  occur  in  this 
meteorological  regime. 

The  definition  of  the  spectral  density  at  time-wave  number  N  [Kahn, 
1957 J  is  the  average  contribution  to  the  variance  from  a  unit  wave  number 
interval  centered  at  wave  number  N.  Since  the  maximum  period  under  con- 
sideration is  50  days,  the  concept  of  a  time-wave  number  N  is  to  be  inter- 


preted  in  terms  of  frequencies  as  follows: 
1  wave  per  50  days  for  N  =  1 , 

N  waves  per  50  days  for  wave  number  N,  N  ^  5. 
Thus  the  frequency  4N    may  be  expressed  as 

fN=JL-  ,      N^  M,  (7) 

and  spans  the  frequency-interval  (N  -  -|-)/2MAt  to  (N  +  i)/2/v1At.  The 
corresponding  period  for  this  time-wave  is  the  reciprocal  of  j-H  .  The 
highest  frequency  ^F  resolvable  when  using  such  discrete  data  at  finite 
intervals  occurs  for  N  =  M  and  is 

j-p  is  often  called  the  folding  or  Nyquist  frequency,  and  in  the  present 
study  is  one  cycle  (1 0  days )  . 

Panofsky  and  Brier  Fl958J  give  the  following  equation  for  raw  spec- 
tral density  estimates  for  time-wave  number  N: 

ttzL 

5N=2/±[R(o)+(-i)MR^ 


1  v  k-1 

In  practice,  SM  is  generally  considered  to  be  distributed  uniformly  over 
the  unit  wave-number  interval  N  -  -|-  to  N  +  •§-,  except  at  N  =  0  and  N  =  M. 
For  the  special  cases  N  =  0  and  N  =  M,  the  corresponding  estimates 
S0  and  Sn  follow  from  equation  (8)  multipled  by  one-half  because  only 
one-half  of  a  wave-number  interval  is  available  at  N  =  0  and  N  =  M. 
These  unsmoothed  raw  spectral  estimates  may  be  smoothed  by  "harming" 
[Blackman  and  Tukey,  1958j,  in  the  following  manner 

S  =  ■  r  S   +  f:  S  4  >  r  S  (9) 
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except  at  the  extreme  wave  numbers  N  =  0  and  N  =  M  where  the  average  is 
modified  to 

3e=.55c+.55i  ,     5m=.5Sm-i+.55m.  (10) 

These  smoothed  estimates  may  be  considered  to  be  centered  at  wave  number 
N,  or  frequency  }n  given  by  equation  (7). 

Note  that  the  integrated  value  of  all  the  spectral  estimates  for  a 
particular  space-wave  n  gives  the  total  time  variance,  for  this  wave.  In 

other  words,  by  definition,  the  relationship 

5 

0ff=£5N>nAN  .    (iD 

holds,  although  carrying  the  integration  only  to  N  =  5  neglects  some  of 
the  higher  frequency  spectral  contributions.  In  addition,  "aliasing"  or 
folding  of  high  frequency  energy  back  upon  the  Nyquist  frequency  has  been 
neglected.  However  the  spectral  densities  at  higher  space-wave  numbers 
(table  2)  suggest  that  this  is  not  important  in  the  present  problem. 
This  has  been  inferred  since  the  contributions  to  the  time  variance  for 
these  wave  numbers  appears  to  act  as  "noise",  with  no  preference  for  the 
higher  time-wave  numbers. 


3.  Computational  procedures  and  interpretations 

As  mentioned  in  section  2,  the  maximum  number  of  time  lags  chosen 
was  M  =  5.  This  choice  gives  a  resolution  of  the  spectrum  of  02./ot 
to  a  maximum  frequency  of  1/10  cycle  (day)  . 

All  computational  processes  were  performed  on  the  CDC-1604.  computer. 
All  programs  devised  for  this  study  were  written  by  the  author  or  modified 
from  those  available  in  the  U.  S.  Naval  Postgraduate  School  computer  cen- 
ter library.  The  three  latitudes  chosen  for  analysis  were  55N,  65N,  and 
75N.  Table  1 ,  a,  b,  c  gives  the  zonally-averaged  autocovariances  for  each 
space-wave  number  computed  by  equation  (6)  at  each  of  the  three  latitude 
circles.  The  smoothed  density  spectra  for  each  space-wave  number  were 
then  computed  according  to  equations  (8),  (9),  and  (10)  and  the  results 
displayed  in  table  2,  a,  b,  c. 

Note  that  in  table  2,  a,  b,  c  periodicities  are  presented  down  the 
left-hand  column,  while  the  space-wave  numbers  n  =  1,  ...,  9  appear  along 
the  top  of  the  table.  From  the  rows  denoted  "P.  R.  Var.",  it  may  be  seen 
that  the  space  waves  n  =  1  and  n  =  2  together  account  for  approximately 
90$  of  the  variance  of  the  zonally-averaged  five-day  height  changes. 
While  the  percent  reduction  of  the  variance  of  dz./d"t  accounted  for  by 
space  wave  2  exceeds  that  of  wave  1 ,  the  contribution  of  the  latter  in- 
creases as  one  proceeds  from  latitude  55^  to  75N.  The  contribution  by 
waves  3  and  J+   accounts  for  less  than  10$  at  latitudes  55N  and  65N,  and 
for  less  than  5%   at  latitude  75N.  Thus,  in  this  study,  only  the  time- 
spectra  of  space  waves  one  and  two  have  been  analyzed  for  significant 
periodicities. 

For  space  waves  n  =  1  and  n  =  2,  fiducial  confidence  limits  at  the 
1 0/o   and  90£  significance  level  have  been  taken  from  Blackman  and  Tukey 
[19531.  These  limits  are  based  upon  the  %/f     distribution,  where 
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(a)  55N 

LAG 
DAYS 

1 

2 

SPACE 
3 

WAVE  NUMBERS 
4             5 

6 

7 

8     . 

9 

0 

42.18 

218.66 

11.30 

13.43 

5.35 

1.64 

2.01 

2.04 

.74 

5 

8.36 

-31 .97 

2.63 

-1.28 

.31 

.14 

.01 

.14 

.04 

10 

2.43 

-31.53 

4.22 

.75 

.26 

-.02 

.16 

.19 

.04 

15 

-3.37 

6.46 

1.71 

2.55 

.79 

-.06 

.32 

.45 

.04 

20 

.86 

5.03 

3.25 

.44 

-.72 

.14 

.40 

.37 

-.03 

25 

-8.46 

-1 1 .47 

2.22 

• 

1.29 
(b)  65N 

.36 

.30 

.32 

.54 

.14 

0 

80.78 

336.38 

15.31 

27.36 

3.11 

.97 

1.37 

1.14 

.27 

5 

23.62 

-139.61 

1.25 

5.23 

.53 

.12 

-.01 

.07 

.01 

10 

-1.55 

25.27 

-2.33 

8.17 

.07 

-.02 

.58 

.09 

-.02 

15 

-4.05 

14.77 

.50 

5.78 

.04 

-.01 

.02 

.18 

-.02 

20 

1.33 

-13.28 

.58 

-.28 

.12 

-.03 

.23 

.24 

.07 

25 

-17.53 

-11.36 

1.18 

3.94 
(c)  75N 

-.10 

.00 

.08 

.11 

.03 

0 

103.59 

139.47 

6.21 

2.69 

1.39 

1.07 

.34 

.36 

.23 

5 

28.08 

-31 .36 

.97 

.44 

-.04 

.22 

.02 

.10 

.04 

10 

7.89 

4.73 

-1.02 

.17 

-.01 

-.03 

-.00 

.03 

.06 

15 

-11.79 

-1.13 

-.10 

.79 

.06 

.09 

.03 

.02 

.02 

20 

-4.05 

42.43 

.19 

-.23 

.25 

.05 

.11 

.04 

-.03 

25 

-11.88 

-46.88 

.41 

-.26 

.01 

.20 

.09 

.04 

-.00 

Table  1.  Zonally  averaged  autocovariances  of  oZ/d't     for  lags  of  0,  5, 
...,  25  days  at  latitudes  (a)  55N;  (b)  65N;  (c)  75N.  The  units 
are  in  terms  of  104  (feet/5  days)2. 


(a)  55N 

FREQUENCY  SPACE  WAVE  NUMBERS 

CYCLES/DAY 

1  23456789 

0/50              9.28     19.34     2.89     1.81        .98       .25       .30  .34  .13 

1/50             10.06     28.31     2.36     1.76     1.00       .31        .29  .31  .14 

2/50              9.30    46.17    1 .61     2.11       .95       .36       .34  .34  .14 

3/50              6.57    54.77     1.40    2.96     1.06       .31       .41  .41  .14 

4/50              4.71     44.96     1.51     2.92     1.01        .24       .36  .36  .13 

5/50 4.05     34.94     1 .61      2.41        .79        .21        .30  .29  .12 

P.  R.  Var.       14.2       77.0      3.5       4.5       1.9       0.7       0.6  0.7  0.3 

(b)  65N 

0/50            18.69     19.95     2.17    6.91        .62       .13       .28  .19  .03 

1/50             20.53     25.75     2.79     5.73       .67       .20       .23  .18  .04 

2/50            19.59     44.51     3.65     3.78       .66       .21        .15  .20  .06 

3/50             13.41     79.03     3.48     3.90       .55       .18       .16  .22  .06 

4/50              7.44     95.69     2.33     4.52       .41        .14       .25  .19  .,04 

5/50 4.38     90.79     1.60     4.13       .33        .11        .29  .15  .03 

P.  R.  Var.       16.0      75.0      3.1       5.2       0.6       0.2       0.2  0.2  0.1 

(c)  75N 

0/50            24.56     13.15       .98       .60       .18       .21        .05  .08  .06 

1/50             26.90     H.79     1.27       .53       .21        .22       .05  .08  .05 

2/50             23.46     24.01     1.58       .46       .27       .24       .07  .08  .03 

3/50            14.54    30.77     1.34       .54       .29       .21        .07  .06  .03 

4/50            10.50     27.46       .81        .46       .24       .14       .05  .04  .04 

5/50 9.57     24.97       .53        .29       .20       .10       .04  .03  .03 

P.  R.  Var.       42.0       53.0       2.6       1.1        0.5          .4       0.1  0.1  0.1 

Table  2.     Zonally-averaged  spectral  densities  for  space-wave  number  n=1 , 
...,9,   and  wave  frequency  N=0,l/50, . ..,5/50  cycles/day.       Sec- 
tions   (a),    (b),   and    (c)  are  for  latitudes  55N,  65N,  and  75N. 
The  units  are  in  terms  of  104(feet/5  davs)2. 
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j—  [Z  '  /  ( —  r1/2 }/  r)   is  the  number  of  degrees  of  freedom.  In  this 
vestigation  we  have  7(  =  35  and  M=  5,  so  that  f  £^£  14- •  For  these  values, 
the  10£  and  90;"  fiducial  limits  interpolated  from  those  given  by  Blackman 
and  Tukey  [l952,  p.20g]  are  0.55  and  1.51  respectively. 

Figure  1  illustrates  the  time-spectral  analysis  at  65N,  for  space 
waves  n  =  1  and  n  =  2  (the  latter  shown  by  the  solid  line).  From  this 
diagram  it  may  be  seen  that  for  n  =  2,  the  mean  variance  averaged  over  all 
periods  is  60.1  x  10  (ft/5  days)  •  The  fiducial  limits  have  been  con- 
structed relative  to  this  mean.  Thus  it  is  evident  that  for  n  =  2,  waves 
of  period  10  days  and  12.5  days  give  significantly  large  contributions  at 
a  90;t  confidence  level.  On  the  other  hand,  at  65N  space  wave  2  at  a  per- 
iod of  50  days,  gives  significantly  small  contributions  to  the  variance. 

The  same  type  of  analysis  for  n  =  1  at  latitude  65N  shows  that  the 
mean  variance  for  space  wave  n  =  1,  averaged  over  all  periods  is  13*5,  in 
units  of  10  (ft/5  days)  ,  so  that  a  time-wave  of  period  50  days  gives  a 
significant  contribution  at  a  90%   confidence  level.  It  should  be  recalled 
that  a  period  of  50  days  in  figure  1  actually  spans  a  range  of  periods 
from  33*3  to  100  days,  and  such  waves  could  include  those  associated  with 
explosive  warming.  This  is  so,  since  strictly  speaking,  the  largest  value 
of  oC,/o X  which  contributes  to  the  variance  of  oZ/ot  could  occur  at 
a  time  which  is  either  three-fourths  or  one  and  three-fourths  of  a  period 
after  the  onset  of  the  autumnal  cooling  season  in  polar  latitudes. 

In  table  3  the  mean  spectral  densities  at  latitudes  55N  and  75N, 
averaged  over  all  frequencies  in  table  2,  are  given  for  space  waves  n  =  1 
and  n  =  2. 


11 


UJ 


< 
»_in 

Q 

uJ 

r 

o 

Q 


s 


< 

ot 

I 


CD 


O 
CD 


O 
00 


t 


O 


■t f 


u 

o 

y- 

Cm 

< 

CO 

Q 

•H 

\ 

VO 

o 

iu 

•H 

^1 

iH 

13 
•H 
H 
O 

<t 

co 

v    Q 

CD 

>■  >_ 

Xj 

o 

E-" 

2 

uj  Q 

• 

3  O 

VA 

LJ  'J- 

CfcUJ 

P 

a) 

LuQl 

• 

*r- 

oJ 

^    II 

o 

(0 

^2 

cd 
> 

in  cd 

\  £ 

•P 

CD    (D 

R»0 

Cm    cd 

^        ' 

w   & 

\CM 

"*       10 

o 

^^ 

v-     M 

o 

Ch  «h 

o 

co 

CO   «H 

°  ^ 

P 

•H    CD 

co  *-< 

CI 

•h  t3 

cd 

CO     r£ 

CD     CO 

•H     Cd 

P  rd 

Q 

•H 

CO    CD 

1  den 
le  th 

<a  .h 

U  A 

p  £ 

o 

CD      »> 

Q<Oi 

o 

CO 

y^o 

II 

CD 

-P    CD 

> 

n3   cd 

X! 

P    CD 

O    O 

o  cd 

tf!    CO 

fO 


C\l  ~ 


o 


O 


12 


LAT 

SPACE  -  WAVE 

1 

2 

55 
75 

7.4 
18.5 

40.4 

23.2 

Table  3.  Mean  spectral  densities  in  units  of  10  (feet/5  days)  at  lati- 
tudes 55N  and  75N  for  space-waves  n  =  1  and  n  =  2. 

Using  the  fiducial  confidence  limits  for  the  X/$     distribution  pre- 
viously mentioned  (0.55  and  1.51,  respectively),  it  is  apparent  from  table 
2  (a,  c)  that  there  exists  no  significant  preferred  peak  in  the  spectral 
densities  at  the  90$  limit  for  either  latitude  55N  or  75N. 

For  space  wave  n  =  1 ,  at  each  of  the  latitude  circles  of  table  2, 
there  is  a  steady  decrease  of  spectral  density,  SN,  with  increasing  time- 
wave  N  5fc  1  • 

For  the  case  of  isotropic  turbulence,  MacCready  [1953J  has  tested  the 
applicability  of  a  spectral  function  of  form 


5NoCN 


-P 


where  N  is  the  time-wave  number.  This  test  is  essentially  based  upon  a 
theory  of  isotropic  turbulence  set  forth  by  Heisenberg  [1948J  in  which 
p  =  5/3.  For  n  =  1,  and  using  the  results  of  table  2,  computation  of  the 
values  of  the  exponent  p  gave  the  following: 

55N,   p  =  0.57 

65N,   p  =  0.91 

75N,   p  =  0.64 
Chiu  n960j  has  obtained  somewhat  similar  results  for  the  exponent  p  using 
the  spectrum  of  the  zonal  and  meridional  components  of  the  wind  at  Belmar 
and  Cocoa,  and  finds  values  of  p  of  the  order  of  1 . 

On  the  other  hand,  the  results  of  this  study  for  space  wave  n  =  2  do 
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not  show  the  decrease  of  spectral  density  with  increasing  time-wave  num- 
ber. Kence  for  n  =  2,  the  peak  spectral  density  which  occurs  in  each  case 
of  table  2  at  either  time-wave  3  or  4-  is  even  more  likely  to  be  signifi- 
cant of  the  average  state  of  synoptic  events. 

Space -wave  n  =  1 ,  with  its  maximum  spectral  contribution  at  time-wave 
N  =  1  is  capable  of  explaining  those  perturbations  of  longer  period  asso- 
ciated with  "explosive"  warming. 
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A*     Conclusions 

(1 )  Space  waves  1  and  2  contribute  9C#  of  the  variance  of  the  height- 
tendency  field  at  latitudes  55N,  65N,  and  75N. 

(2)  The  chief  contributory  to  the  spectrum  of  oZ/ot   occur  for 
time-graves N  =  3  or  N  -   4.  (periods  12.5  and  16.7  days),  which  have  the  same 
general  periodicity  as  the  average  perturbation  period  at  a  fixed  station 
[see,  for  example,  Austin  and  Krawits,  19561. 

(3 )  The  chief  contributor  to  the  time  variance  for  space  wave  n  =  1 
occurs  for  time  wave  II  =  1 ,  or  possibly  l\  -  -tr,   which  corresponds  to  periods 
of  the  order  of  50  days  or  longer.  Such  periods  may  be  sufficiently  long 
to  account  for  the  "explosive"  warming  which  sometimes  occurs  as  early  as 
January. 

(4.)  For  space  wave  n  =  1,  the  time  spectral  density  was  distributed 
with  time  wave  number  according  to 

r    k  i  -a?l 

SNc£N 

which  compares  reasonably  with  the  results  obtained  in  other  spectral 
studies  fc.f.  Chiu,  196ol. 
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APPENDIX  I 
DERIVATION  0?  THE  ZONALLY-AVERAGED  TIME  AUTOCOVARIANCE 
The  equa.tion  for  the  zonally-averaged  tine  autocova.ria.nce  of  the 
height-tendency  field  at  lag  X    may  be  expressed  in  the  form  shown  in 
equation  GO 


^_  ,  -    —   x  ^   _  ,j   _ .  -  7 

where  Tj  is  the  maximum  nunber  of  discrete  values  in  the  sample,  and  R(k) 
is  the  zonal ly-averaged  time  autocovariance . 

When  equation  (4. )  is  expanded  using  equations  (2)  and  (5),  and  upon 
interchanging  the  order  of  summation,  equation  (12)  is  obtained 

^M         /2TT 

Rn(K)  ^V±  |    (Mi[cos(2nA-(()n-(pn)  +  COo(M^)]+ 


Cn4i  [sin  (2  n  A  -  (j)h  -%)  +  sin  ( l|*r  0n )] 


+ 


q^Cn  O  f 2 nA  -  '\  -  ip )  - s-M  ip,:  -0,)]  +       I 

.    &  I,  W,     (C0s(lp,   -0ni-COs(2nA-([)n-^.)]ji/  J,(12) 

where  the  orthogonality  property  has  also  been  used.  Equation  (12)  repre- 
sents a  double  integral;  first,  with  respect  to  longitude,  second,  with 
respect  to  time.  The  integral  with  respect  to  longitude  gives  the  zonal 
average.  On  integrating  around  the  latitude  circle  from  0  to  27T,  all 
cosine  and  sine  terms  involving  the  argument  (2  HA  —  (£>n — (J)A] 
integrate  to  zero.  Therefore,  the  result  obtained  is 
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CJ 


'/.  K 


+ 


R(K)-)  _l  ;   ,  c,i.  '-c,,g^\^,]r-  fiji,-^ ) 

fo'lH^l TT 


(6) 


i'ote  that  nn(K)   is  the  sonally-averaged  antocovariance  at  a  time-lag 
of  K&t  days,  fcr  the  nth  space-wave;  no  other  space-wave  appears  in 
the  expression  for  F?n(K)  •  'The  interval  between  successive  data  is 
such  that  A  L  —  5    days . 
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